Scaling and memory in the return intervals of energy dissipation rate in 
three-dimensional fully developed turbulence 
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We study the statistical properties ol return intervals r between successive energy dissipation 
rates above a certain threshold Q in three-dimensional fully developed turbulence. We find that 
the distribution function Pq(t) scales with the mean return interval Rq as Pci(r) = Rq 1 f(r/Rg) 
except for r = 1, where the scaling function f(x) has two power-law regimes. The return intervals 
are short-term and long-term correlated and possess multifractal nature. The Hurst index of the 
return intervals decays exponentially against Rq, predicting that rare extreme events with Rq — > do 
are also long-term correlated with the Hurst index Hoc = 0.639. 
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Extreme events are ubiquitous in nature and society 
and understanding their dynamics is of crucial impor- 
tance 0, However, extreme events are usually rare, 
which makes it difficult to investigate their occurrence 
properties. Recently, there is an increasing interest in 
the study of return intervals or reoccurrence times r be- 
tween successive events above (or below) some threshold 
Q, aiming at unveiling the laws governing the occurrence 
of extreme events by studying the statistics of the return 
intervals for increasing threshold Q. 

Recent studies show that the long-term correlation 
structure has essential influence on the statistics of re- 
turn intervals fsj]. For long-term power-law correlated 
monofractal records with exponent 7, numerical anal- 
ysis illustrates that the distribution density function 
of the return intervals follows a stretched exponential 
Pq(t) ~ exp [— b(r / Rq) 1 ] with the same exponent 7 and 
the return intervals are also long-term correlated, again 
with the same exponent 7, where Rq is the mean re- 
turn interval associated with threshold 

Q [Illll- Un- 
corrected records with 7 = 1 is a special case, whose 
return intervals arc exponentially distributed. When 
< 7 1, theoretical analysis with certain approxi- 
mation shows that the distribution of return interval is 
Weibull: P Q (r) ~ r^ 1 exp(-crT) This seems con- 

sistent with recent numerical results [9( . 

For multifractal records in the presence or absence of 
linear correlations, extensive simulations based on the 
multiplicative random cascade (MRC) model |10| and 
the multifractal random walk (MRW) model [11| un- 
veil that the return intervals have a power-law decay in 
the distribution and are long-term correlated governed 
by power laws whose exponents depend explicitly on the 
threshold Q, and the conditional return intervals increase 



as a power-law function of the previous return interval 
13, 0, Q. These results are of particular interest since 
a variety of time series exhibit multifractal nature. For 
instance, the returns of common stocks can be well mod- 
eled by the multifractal random walk 15, [HI], and the 
statistics of the associated return intervals are found to 
comply with the numerical prediction l^, | , which can 
be used to significantly improve risk estimation [l2l ]. 

Another important issue of multifractal records con- 
cerns the possible scaling behavior of the return in- 
terval distributions -Pg(r) over different thresholds Q. 
Numerical simulations of MRC and MRW time series 

However, 
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find no evidence of such scaling [12l . |R 
empirical return interval analysis of financial volatility 
gives miscellaneous results. Several studies reported that 
there is a scaling law in the return interval distributions 
17 , H, H, 2(| 21 1 , while others argued that the cumula- 
tive distributions of return intervals had systematic de- 
viations from scaling and showed multiscaling behaviors 

H SHI Hi- 



In this Letter, we perform return interval analysis of 
the energy dissipation rate in three-dimensional fully de- 
veloped turbulence based on a high-Reynolds turbulence 
data set collected at the SI ONERA wind tunnel by the 
Grenoble group from LEGI [26j. The size of the veloc- 
ity time series {vi : i = 1, 2, ■ • ■ , N} is about 1.73 x 10 7 . 
Using Taylor's frozen flow hypothesis which replaces a 
spatial variation of the fluid velocity by a temporal vari- 
ation measured at a fixed location, the rate of kinetic 
energy dissipation at position i is e; ~ [(«;+i — Vi) /5e] 2 , 
where 81 is the resolution (translated in spatial scale) of 
the measurements. The energy dissipation rate time se- 
ries exhibit multifractal nature (27l | and its Hurst index 
is H = 1 — 7/2 = 0.81. Contrary to previous studies, 
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we find that the return interval distributions show two 
power-law regimes and collapse onto a single curve for 
different thresholds Q. The scaling phenomenon is also 
observed for conditional interval distributions. 

We have calculated the return interval time series r, for 
different thresholds Q, which can be mapped nonlinear ly 
to the mean return intervals Rq . Logarithmic binning is 
adopted to construct the distribution density functions 
Pq(t). In order to ensure that the bins cover the whole 
r-axis, we use the following procedure. First, the inter- 
val [l,max(rj)] is partitioned logarithmically into n — 1 
subintcrvals whose edges are X\ < x 2 < ■ • • < x n . Then 
we obtain the sequence yi = [xi\, where [ Xj^ IS Si round 
function of Xi. We discard duplicate integers in the yi se- 
quence and obtain a new sequence Wj . The edge sequence 
of the bins are determined by {&;} = {0.5, {wi + 0.5}}. 
For each bin (e^, e^+i], the empirical density function can 
be calculated by 



P Q (n) = 



#(e 4 < r < e i+ i) 



f 



#{r > 0) e l+1 



(1) 



where = (e^ + ej+i)/2 and #() is the number of return 
intervals that satisfies the condition in the parenthesis. 
The empirical distribution of the return intervals is de- 
picted in Fig. Q] for Rq = 50, 150 and 500. We find that 
the three distributions collapse onto a single curve ex- 
cept for r = 1, in remarkable contrast to the simulation 
results for MRC and MRW time series. Two power-law 
regimes are observed 



RqPq(t) 
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where A\ = 0.107 and A 2 = 33.4 are prefactors, the 
crossover return interval r c m 7Rq, and S\ — 0.987 ± 
0.013 and S 2 = 3.88 ± 0.09. For the shuffled data, the 
RqPq(t) curves collapse to a single exponential curve. 

In risk estimation, a quantity of great interest is the 
probability Wc}(At, t) that an extreme event occurs af- 
ter a short time At <C t from now on, conditioned that 
the time elapsed t after the occurrence of the previous 
extreme event 
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j t t+At p Q m 
f t °°p Q (t)dt ■ 



W Q {At\t) 



(3) 



When t > r c . simple algebraic manipulation leads to 



W Q (At\t) ~ (82 - l)At/t. 



(4) 



The probability Wq(At\t) is found to be proportional to 
At and inversely proportional to t. An intriguing fea- 
ture is that Wq(At\t) is independent of the threshold Q, 
which is a direct consequence of the scaling behavior of 
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FIG. 1: (color online). Log- log plot of the scaled distribution 
RqPq(t) as a function of t/Rq for three different values of 
-Rq. The straight lines are the best fits to power laws. The 
power-law exponents are S\ — 0.987 ± 0.013 and 82 = 3.88 ± 
0.09. For the shuffled data, the RqPq{v) curves collapse to a 
same exponential curve, as expected. 



Pq(t) shown in Fig. [TJ When t < r c , we obtain that 
W Q (At\t) 
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We find that Wc}(At\t) is proportional to At. However, 
Wq(At\t) also depends on Rq. 

In order to test the memory effects of the return inter- 
vals, we first investigate the conditional PDF Pg(r|ro), 
which is the distribution of return intervals immediately 
after r . To gain better statistics, we study PQ(r\ro) 
for a range of ro rather than individual ro values. For 
each threshold Q or Rq, the return intervals sequence 
are sorted in an increasing order and then divided into 
eight groups G±, ■ ■ ■ ,Gg with approximately equal size. 
An empirical conditional distribution is determined for 
each r group. Figure shows Pg(r|r ) for r G G\ and 
ro eGg. For each group, the three distributions collapse 
onto a single curve, indicating evident scaling behavior. 
The figure shows that the probability of finding small 
(large) r in G\ is enhanced (decreased) compared with 
Gg . This discrepancy in the two groups of distributions 
unveils the memory effect that large (small) return inter- 
vals tend to follow large (small) return intervals. This is 
true since the distributions associated with different Gi 
should not exhibit significant discrepancy if there is no 
memory in the return intervals 17 1. 

The memory effect in the conditional distribution 
PQ(r\ro) can also be illustrated by the mean conditional 
return interval (r\ro). If there is no memory in the return 
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FIG. 2: (color online). Scaling and memory in the conditional 
distributions. The scaled conditional distribution R.QPQ{r\ro) 
is plotted as a function of r /Rq with ro in Gi and Gs for three 
different values of Rq. 



intervals, (r|ro) does not depend on ro and is equal to Rq. 
Figure [3] plots the mean conditional return interval (r|ro) 
as a function of ro for Rq = 50, 150 and 500. It is shown 
that (r|ro) is a power-law function of ro when ro is larger 
than certain value and the power-law exponent decreases 
with Rq , which is consistent with many simulational and 
empirical studies. Also shown in Fig. [3] is the mean con- 
ditional return interval of the shuffled energy dissipation 
rate, which does not depend on ro. 
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FIG. 3: (color online). Scaled mean conditional return in- 
terval {r\ro)/RQ versus ro/Rg of return interval sequence in 
turbulence for various Rq. The filled symbols represent the 
results for the shuffled records without memory. 



mial function and the local rms fluctuation fk(s) is de- 
termined. The overall detrended fluctuation is calculated 
by 



1/9 



(0) 



By varying the value of s, one can expect the detrended 
fluctuation function F q (s) scales with the size s: 



F g (s) 



Mq) 



(7) 



where h(q) is the generalized Hurst index. The return 
interval series possesses multifractal nature if and only 
if h(q) is a nonlinear function of q. When q = 2, h(2) 
is nothing but the Hurst index H and the MFDFA re- 
duces to the DFA. The Hurst index H is related to the 
autocorrelation exponent 7 by 7 = 2 — 2H. 

We first investigate the linear long-term correlation 
property of the return intervals using DFA. The depen- 
dence of the fluctuation function F% is drawn in Fig. [4] 
against s for Rq = 50, 150, and 500. In all cases, we find 
nice power-law relation and the scaling range decreases 
with the increase of Rq. The inset shows the depen- 
dence of the Hurst index with respect to Rq, which has 
an exponential decay: 

H = H 00 + be- RQ / R ° = 0.639 + 0.158e~ R » /69 - 9 , (8) 

where R c = 69.9 is the characteristic scale. For extreme 
events with very large Q, Rq tends to infinity, and the 
Hurst index can be predicted as H = = 0.639. This 
implies that the return intervals of those extreme events 
also exhibit long-term memory. 
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We now study the long-term correlation in the re- 
turn intervals using the multifractal detrended fluctua- 
tion analysis (MFDFA), which is able to extract long- 
term power-law correlation in non-stationary time series 
HSHEHfll]. The MFDFA considers the cumulative time 
series Ri = X^=i( r » ~ ( r )); which is partitioned into N s 
disjoint boxes with the same size s. In each box fc, the 
local trend is removed from the subseries by a polyno- 



FIG. 4: (color online). Detrended fluctuation analysis of the 
return interval time series for different Rq. The inset shows 
the exponential decay of the Hurst index H against Rq. 

We also apply the MFDFA to investigate the multi- 
fractal nature of the return intervals. Figure [5]^a) shows 
the power-law dependence of the overall fluctuation F q (s) 
on the scale s for Rq — 50. When Rq increases, the 
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F q (s) function becomes more noisy, especially for neg- 
ative q. The slopes of the straight lines are the linear 
least-squares estimates of the generalized Hurst indexes 
h(q), which are drawn in Fig. [5jb) . The mass scaling 
exponent function r(q) = qh(q) — 1 and the multifrac- 
tal spectrum f{a) calculated according to the Legendre 
transform of r(q) are illustrated respectively in Fig. [5tc) 
and Fig. (5^d). The sound nonlinearity in h(q) and r(q) 
is a hallmark of multifractality in the return intervals, 
whose singularity strength increases with Rq. 
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FIG. 5: (color online). Multifractal detrended fluctuation 
analysis of the return interval time series for different Rq. 
(a) MFDFA fluctuation function for different q when Rq = 
50. (b) Generalized Hurst indexes h(q). (c) Mass scaling 
exponents r(g). (d) Multifractal singularity spectra /(a). 

In summary, we have studied the statistical proper- 
ties of return intervals of the energy dissipation rate in 
three-dimensional fully developed turbulence. The in- 
terval distribution is found to exhibit scaling behavior 
across different Rq and two power-law regimes, except 
for intervals r = 1. We found that the conditional in- 
terval distributions also collapse onto a single curve for 
same tq, but deviate for different Tq, and the mean con- 
ditional interval increases as a power law of ro, which 
indicates the presence of short memory in the return in- 
tervals. The long-term memory and the multifractal na- 
ture in the return intervals are confirmed by the DFA 
and MFDFA. The Hurst index of the return intervals de- 
cays exponentially against Rq , which allows us to predict 
the asymptotic Hurst index of return intervals between 
rare extreme events as Hoc = 0.639. Our results signify 
discrepancy with the numerical results using MRC and 
MRW models [13], implying that the energy dissipation 
process can not be modeled using these models. 
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